Dynamics of Spin Relaxation in Finite-Size 2D Systems: an Exact Solution 
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We find an exact solution for the problem of electron spin relaxation in a 2D circle with Rashba 
spin-orbit interaction. Our analysis shows that the spin relaxation in finite-size regions involves 
three stages and is described by multiple spin relaxation times. It is important that the longest 
spin relaxation time increases with decrease in system radius but always remains finite. Therefore, 
at long times, the spin polarization in small 2D systems decays exponentially with a size-dependent 
rate. This prediction is supported by results of Monte Carlo simulations. 
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The problem of D'yakonov-Perel' [HE] spin relaxation 
in two-dimensional (2D) systems have attracted wide at- 
tention flT-TH] because of its fundamental importance for 
the field of spintronics [151 [TB] , However, the spin re- 
laxation in systems with boundaries is even more im- 
portant because boundaries are naturally present in all 
electronic devices. There are only several examples in 
the literature where the influence of boundary conditions 
on D'yakonov-Perel' spin relaxation have been explored 
theoretically and/or experimentally. These examples in- 
clude investigations of spin relaxation in 2D channels |17f - 
[22] . 2D half-space [23], 2D systems with antidots [24], 
large quantum dots |25l |2"5], and one-dimensional (ID) 
finite-length wires [57]. Both available experimental and 
theoretical results indicate that typically in the diffusive 
spin transport regime the electron spin life time is longer 
in systems with boundaries. 

In this paper we find an exact solution for the prob- 
lem of electron spin relaxation in finite-size 2D systems. 
Specifically, we consider the dynamics of electron spin 
relaxation in a 2D circle made of a semiconductor struc- 
ture with Rashba-type spin-orbit interaction. Naively, 
one may think that in small systems the spin relaxation 
is incomplete as the spin precession angle across the sys- 
tem is small. However, in such a situation, the different 
effect plays role: the non-commutativity of spin rota- 
tions. Because of this effect, the electron spin precession 
angle can largely exceed the maximum rotation angle al- 
lowed by naive geometrical considerations. To the best 
of our knowledge, the spin relaxation in small systems 
was investigated previously only in Ref. 1251 This previ- 
ous study provides only an asymptotic value of the spin 
relaxation time without giving details on how the whole 
process of spin relaxation occurs. In the present paper, 
we show that the spin relaxation process in finite-size sys- 
tems is intrinsically complex. The exact solution of this 
fundamental problem involves an infinite number of spin 
relaxation constants. At long times, however, only the 
slowest decaying component survives and the spin po- 
larization exhibit a slow size-dependent exponential de- 
cay. Our exact analytical solution is obtained using the 



Laplace transform and is confirmed by Monte Carlo sim- 
ulations of spin dynamics. This work thus provides an 
important missing part of spin relaxation theory. 

Let us consider the dynamics of electron spin polariza- 
tion in a finite-size 2D electron system such as the circle 
shown in Fig. [T] The set of diffusion equations [H [23 121] 
for spin polarization is given by 

= AS in + 2 V WS Z - r? 2 S in , (1) 

1 dS 

-^=AS Z -2 V V-S in -2 V 2 S Z , (2) 

where Si n and S z are the in-plane and z components of 
spin polarization, respectively, D = 1 2 /(2t) is the diffu- 
sion coefficient, I is the mean free path, r is the momen- 
tum relaxation time, and rj is the spin rotation angle per 
unit length (spin rotations are induced by the Rashba 
spin-orbit interaction [30] )• The applicability limits of 
Eqs. jl} and (|2|) are A <C I <C a and Z <§C 77 1 , where 
A is the electron's de Broglie wavelength, and a is the 
characteristic system size (in the context of this work, a 
is the radius of the circle). We also note that Eqs. Q 
and Q are valid for any value of arj. Eqs. ([IJ and Q are 
supplemented by the standard boundary conditions |28| 
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FIG. 1: (Color online) Schematic of a circular 2D region of a 
radius a with electron spin polarization initially pointing in z 
direction perpendicular to the circle plane. 
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where n and T are the normal in-plane and tangential 
vectors to the boundary T, respectively. 

Next, we would like to reduce the set of Eqs. (JTJ and 
([2]) to a single equation for S z . Introducing u = V • S; n 
and v = (V x Sj n ) z , Eqs. ([!]) and Q can be rewritten as 



u t = Au + 2AS Z — u, 
(S z ) t = AS Z - 2S Z - 2u, 
v t = Av — v. 



(4) 
(5) 
(6) 



Here, A is the 2D Laplace operator, the time is measured 
in the units of t s — (Dtj 2 ) -1 , and the coordinates are 
measured in the units of Such a convention is used 
below if not stated otherwise. Combining Eqs. Q and 
([5]) we readily get 

(S z ) tt - 2A(S z ) t + 3(S z ) t + A 2 S Z + AS Z + 2S Z = 0. (7) 

The Laplace transform of Eq. ^ is given by 

A 2 S Z + (1 - 2p)AS z + ( P 2 + 3p + 2)S Z 

= (p + l)S z \t=o-AS z \t=o-2u\ t = , (8) 

where S z is the Laplace transform of S z to the complex 
p-domain. In the above equation, the time derivative of 
S z at t — was substituted from Eq. ([5|. 

In what follows we consider the relaxation of homoge- 
neous spin polarization pointing in z direction perpen- 
dicular to the plane of 2D circle (see Fig. [IJ. Taking into 
account the axial symmetry of the problem, the initial 
and boundary conditions read 



u(r, t = 0) = 

dS r 
dr 



d(rS r (r,t = 0)) 



rdr 



= 0, 



--0,S z (r,t = 0) = S , (9) 
dS z . 



= 0. (10) 



We note that the function v = d(rS^)/(rdr) is safely 
taken out of the consideration since the solution v(r, t) = 
satisfies Eq. (|6| with the boundary condition 
dScf,/ 'dr\ r=a — and the initial condition v(r, 0) = 0. 
Applying the initial conditions Eq. Q simplifies to 

A 2 S z (r) + (l-2p)AS z (r) + (p 2 +3p+2)S z (r) - (p+l)S . 

(11) 

The general solution of Eq. (Ill can be found by the 



factorization of its left-hand side and be presented as 



S z (r) = A 1 J {k 1 r) + A 2 J (k 2 r) + BMfar) 

So 



+B 2 N (k 2 r) + 



P 



2' 



(12) 



where Ai^ 2 ,B\ t2 are arbitrary constants, Jq(z) and Nq(z) 
are the zeroth order Bessel and Neumann functions, re- 
spectively, and 
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(13) 



In the case of the circle, B\ 2 = since Nq(x) diverges 
as x — > 0. Moreover, the actual choice of two branches 



corresponding to ± in Eq. (13) is not essential. We may 



make a branch cut along the line p < —7/16, Im(p) = 
in the plane of complex p and define two branches by the 
conditions k\ 2 {p = 0) = (1 ± iy/7)/2. 

Although the radial component of spin polarization 
equals zero at t = 0, it becomes different than zero at 
t > similarly to the case of spin relaxation in rings |31j . 
With a help of Laplace transform of Eq. (|5| , we express 
S r (r) through S z (r) as 



2r 



S r (r) = - 
1 

d^SAO 



ies z (t 



1 
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2 dr 
dttS z (£,t = 0)- (14) 



The boundary conditions ( 10 ) are used to find the val 



ues of Ax 2 in Eq. (12). For this purpose, we Laplace 



transform Eqs. (10 I and employ Eq. (14) to find 



d 2 S z (a) p + 2 
dr 2 a 2 



drrS z (r) - pS x (a) 



1 f a 

T-\ I drrS z (r,t = 0) - SJa,t 
a 2 Jo 



0),(15) 



AS, (a) | p + 2 
dr a 



drrS z (7 



drrS z (r,t = 0).(16) 



10 u Jo 

A\, 2 are obtained from Eqs. (15) and ( |16[) c omplemented 
by Eq. (12 1 and the initial conditions (|9|. Finally, the 
Laplace transform of z-component of spin polarization is 
written as 
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(17) 



where the following notation is used: 
D(p) = 2(p + 2)(k 2 -k 2 ) 



ak\k 2 



+ [(p + 2)(k 2 1 -l)-pk 2 2 }Mk 1 a) 
~[{p + 2){k 2 -l)-pk 2 ]Jo{k 2 a) 



• , 2 ^ Ji(fcia) Ji(fc 2 Q) 
ik 2 

Ji(k 2 a) 

k 2 
Ji(fcia) 



(18) 



The Laplace transform of r-component of spin polariza- 

J\{k 2 a) 



tion is found combining Eqs. (14 1, (17) and (|9 

So 



Sr{r) - (p + 2)D(p) 
Ji(kir) n2 



{ki+p + 2){k 2 2 -p-2) 



(k 2 +p + 2)(ki-p-2) 



k 2 

Ji(fcia) Ji(fc 2 r) 
h k 2 



.(19) 



The inverse Laplace transform of Eqs. (17) and (19 1 



provides the time-domain components of spin polariza- 
tion. It is important to note that the right-hand sides of 
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FIG. 2: (Color online) First six poles of S x (r) and5V(r) versus 
the ring radius. These poles are proportional to the slowest 
six spin relaxation rates. 



Eqs. (17 1 and (19 1 do not change under a permutation 
of k\ and fc 2 . This means that these functions are one- 
valued functions in the wh ole complex plane of p despite 
the square root in Eq. (13 1. As a result, S z (r) and S r (r) 
are meromorphic functions. The poles of these functions 
are defined by the equation D(p) = 0, which have an in- 
finite number of roots p n , \p n \ — > oo as n — > oo. All poles 
are characterized by Im(p„) = and Re(p„) < 0. Note 
that, generally, p = —2 is not a pole of both S z (r) and 
S r (r) . Referring to Fig. [5J the positions of poles depend 
on the circle radius. At small values of arj, the values of 
all \p n \ (n = 1, 2, ...) increase with decrease of arj except 
of \po\ whose value decreases. Basically, this is the most 
interesting pole describing the asymptotic spin relaxation 
at long times, po is always located between —7/16 and 
and tends to zero when arj — > 0, and to —7/16 when 
arj — > oo. In the limit of small a, an analytical expression 
for p can be found. Expanding Eq. (18 1 over small p 
and a, we get (using dimensional units for a) 

1 , ,A 7 



Po 



4cS 
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(ar)) + O ( (arj) 



(20) 



This expression is basically valid when arj < 1. We also 
note that the first term in the right-hand side of Eq. ( 20 ) 
was previously reported in Ref. 1251 It can also be shown 
that for n = 1,2, \p n \ ~ C n (arj)~ 2 when arj <C 1. 

The time-domain components of spin polarization are 
written using the normal dimensional units of time and 
coordinates as 



+°° QO „p n Dn 2 t 

Sz(r,t) = 

n=0 



( Pn + 2)D'(p n 



(Pn 



, 2 ^Ji(k2 n arj) 
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(21) 
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FIG. 3: (Color online) Radial distributions of S r and S z com- 
ponents of spin polarization at different moments of time. 
This plot was obtained at arj — 1. 



and 



+ °° Q pPnDrft 
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ki 



(22) 



Eqs. (21 1 and (22) represent the main result of this work 
describing the time-dependence of spin relaxation in the 
circle. Basically, three main stages of spin relaxation in 
small systems can be identified (this separation is appro- 
priate at arj < 3, when, as it follows from Fig. [2j p$ is 
well separated from all other poles). The first (initial) 
stage of spin relaxation takes place at t < a 2 /(16D), 
where a 2 /(AD) is the time it takes for an electron to dif- 
fusively propagate over a distance a. During this stage, 
the most of electrons in the circle's center still do not 
'know about' the presence of the boundary and, there- 
fore, the spin relaxation occurs essentially as in the bulk 
(accordingly to the standard 2D D'yakonov-Perel' spin 
relaxation theory). In the second stage of spin relaxation, 
when a 2 /(1QD) < t < &/(\pi \Drj 2 ), several exponentially 
decaying terms play the main role in Eqs. (21 ) and (22 1. 



During this stage, a slow-decaying spin polarization pro- 
file establishes. In such a profile, \S r \ increases with r as 
we demonstrate in Fig. [3] The last third stage of spin 
relaxation is a slow single-exponent decay of the slow- 
decaying spin polarization profile established during the 
second stage. This process occurs at long times, namely, 
when 4/(|pi \Drj 2 ) < t. The analytical expressions de- 
scribing the shape of the slow-decaying spin polarization 
profile can be easily inferred from Eqs. (21) and (22 1 in 
the long-time limit. All three stages of spin relaxation 
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FIG. 4: (Color online) Time dependence of S z component of 
spin polarization at r — obtained using Eq. (21 1 and Monte 
Carlo calculations. The dashed line corresponds to the usual 
D'yakonov-PereP relaxation in infinite 2D systems. This plot 
was obtained using the parameter value rjl = 0.1. 



can be easily distinguished in Fig. |4j 

In order to obtain an additional insight on spin relax- 
ation in 2D circle, we have performed extensive Monte 
Carlo simulations. All specific details of the Monte Carlo 
simulations approach can be found in Refs. [17] and [32] 
and will not be repeated here. We just mention that 
the Monte Carlo simulation program uses a semiclassi- 
cal description of electron space motion and quantum- 
mechanical description of spin dynamics. A spin conser- 
vation condition was used for electrons elastically scat- 
tered from system boundaries. Generally, all obtained 
Monte Carlo simulation results are in perfect quantitative 
agreement with our analytical predictions thus confirm- 
ing our analytical theory of spin relaxation in finite-size 
systems. A comparison of selected analytical and numer- 
ical curves is given in Fig. [4] 

In summary, we have found an analytical solution for 
the problem of electron spin relaxation in the circle. It is 
shown that a small but non- vanishing spin relaxation ex- 
ists even in small systems at long times. Consequently, 
it is not possible to completely eliminate the electron 
spin relaxation by reducing the system size, although the 
relaxation rate is dramatically suppressed in small-size 
systems. Basically, the spin relaxation process can be 
separated into three stages including an initial region of 
fast bulk- type relaxation, a transition region where the 
relaxation is described by a combination of several ex- 
ponentially decaying functions and a region of slow ex- 
ponential decay at long times. Our results can be easily 
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